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Abstract 

A large class of conformal field theories can be studied by setting the mass to zero in a set 

of non-interacting fields that form a representation of the Lorentz group. We consider all such 
O 

field theories defined on a general compact manifold. In the large volume limit, the density of 
states as a function of energy is shown to obey the same differential equation regardless of the 

(N 

manifold's geometry. From this equation, we compute the leading and next-to-leading behaviour 
for the density of states, and establish a method that can be iterated to obtain the asymptotic 

Ph 

<j_) series to any desired precision. 
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I. INTRODUCTION 



The density of states is an example of a thermodynamic quantity which, even in a free field 
theory, is only ever calculated approximately. Part of the problem is that the notion of 
a measure p(E)dE giving the number of states with energy in the interval [E, E + dE) is 
itself an approximation. For a theory with a discrete spectrum, this assumption is only valid 
when the energy is high. For this reason, it is natural to approximate the density of states, 
asymptotically. The goal of this work is to better understand the asymptotics of p(E) in a 
free conformal field theory. 

Conformal field theories in 1 + 1 dimensions have been studied extensively since the dis- 
covery of worldsheet string theory. More recently, Maldacena's discovery of the AdS/CFT 
correspondence [1] has ignited an interest in higher dimensional CFTs[2H6] as well. Our re- 
sults hold for a general d where d is the number of spatial dimensions. For d > 2, the results 
appear to be new. 

Much is known about the density of states in a two-dimensional CFT (d = 1 in our notation). 
By convention, the volume in such a theory on a circle is usually normalized to 2ir. A famous 
result is the Cardy formula[7j which expresses the density of states in terms of the central 
charge c (defined so that the vacuum has an energy of — ^). It follows from Cardy's 1986 
paper that 1 

p(E) oc e 27r Vf = 

for large E. Clearly, this cannot be an exact formula. It does not have an inverse energy as 

1 For an explanation of the notation, we will use oc in equations that are roughly true to illustrate a 
conceptual point. For example f(x) oc g(x) means that the leading term of / is proportional to the leading 
term of g. A symbol we use more rigorously is ~. We say that f(x) ~ g(x) if f(x) = g(x) + o(g(x)). In 
other words, the leading terms are equal and as we let x — > oo, / — g becomes negligible compared to / 
and g. 
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a prefactor for instance. In fact, in 2011, Loran, Sheikh-Jabbari and Vincon showed that 



p{E) 



3 2n v /cE/3 



rccV J i ^cVE/Z 



3 y^cVE/3 

where I\ is the modified Bessel function of the first kind[8j. Bessel functions have a well 
known asymptotic series [9] and we can use this series to see that the above formula is 
consistent with and extends the Cardy formula. For example, the crudest approximation 
beyond the Cardy formula is 

p{E) ^ w e^^ (1) 

6(27r)3(d/£/3)! 

The main property assumed in these derivations is that the partition function is invariant 
under modular transformations. For our starting point, we instead assume that the partition 
function is that of a free theory. This allows us to treat the case of higher dimensions in 
which a manifold is not necessarily flat. 

Using dimensional analysis and conformal invar iance, there is a simple argument that reveals 
how p(E) must roughly scale with E for a translation invariant system. Energy is extensive 
so it must be proportional to the volume. For a CFT, the only other scale is the one set by 
the temperature so by dimensional analysis, 

E = CVT d+1 

where C is some dimensionless constant. The same arguments can be applied to entropy, 
giving us: 

S = CVT d ocV^E^ 
Finally, Boltzmann's formula relates the entropy to the density of states: 

d 

p{E) oc e oc e 



The goal of this work is to improve on this formula. If ((a) = Y^=% ^ * s the Riemann zeta 
function, C*(°") = J2™=i — i s ^ ne alternating zeta function and Ud is the volume of a 
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unit ball in R , the result that will be derived is: 



(d+l) d +M d +! 



a d+i e d+i 



p(£) ~ CjaWEWel 7 (2) 

rf!u; d [sC(rf + l) + s*C*(rf+l)]^ 
° ~ (2vr) d 

in a theory that has s species of bosons and s* species of fermions. Cd is a constant depending 
only on the dimension which can be computed without much difficulty When only the 
leading behaviour of the density of states is desired, it is common to see the entropy written 
instead. Taking the logarithm of our main result, 

logp(E) ~ ( (d+ J, )d+1 ) d " a^E^ - log (C d aE^) (3) 

The prefactor in ^ or equivalently the log term in ^ is new and will be referred to as 
"the logarithmic correction" in this work. The assumptions we make about the physics are 
captured in the value of a. Once a is known, the form of ^ is generally true. Consider, for 
example, the maximally supersymmetric large N Yang-Mills theory. This theory is widely 
studied for its direct role in the AdS/CFT correspondence. The Super Yang-Mills entropy 
for a coupling parameter of zero 

n f25Qir 2 N 2 V\ * » 

*=Ur— ) EI 

appears frequently in the literature [10HT3j. Substituting d = 3 and s = s* = 8N 2 into our 
formula for a, we see that the entropy above is precisely the first term in ([3]). We have 
therefore computed a logarithmic correction to the well known result. In the interacting 
version, [TUHT2] have found that the physics can still be well described if the non-interacting 



a value of ir2N ^ v is multiplied by a function of the t'Hooft coupling parameter (which is | 
in the infinite volume limit). Our analysis does not have anything to say about what this 
funciton might be, but once it is specified, terms in the asymptotic series for p(E) can be 
found using the same method that we apply to a free theory. Another non-free theory in 
which our result might hold is quantum chromodynamics in the high energy limit. This non- 
conformal field theory has a scale Aqcd ~ 217MeV which emerges at the quantum level. A 
famous result [TU [15] is asymptotic freedom which states that the running coupling of QCD 
vanishes at the highest energy scales. Specifically, 

1 



where q 2 is the squared momentum and rif is the number of fermions. A study of the 
thermodynamics in this limit [TB] has found that the free energy is given by 

8ttV 



45 



-T 



Their first term is that of a free theory, consistent with asymptotic freedom. If we use only 
this first term of the free energy, we can derive 



b ~3 



ln 2 V ( 21 

1 + ™ n f 



15 V 32 

which is the first term in our ^ with d = 3, s = 16 and s* = 12nf. However, there are 
corrections to the free energy coming from the fact that QCD is not free. These corrections 
that depend on q 2 (and therefore on the temperature) cannot simply be absorbed into the 
value of a. Unlike with Super Yang-Mills, there is no simple adjustment to ^ that would 
allow it to account for the coupling constant. One would have to already know that it is 
close to zero a priori. 

The derivation will be done in three main steps. In the first step, our thermal partition 
function is calculated on M x {time} where M is a compact <i-dimensional manifold. This 
manifold may or may not have a boundary. In the second step, the density of states is 
expressed as a Taylor series which is found by taking the inverse Laplace transform of the 
partition function (transitioning from a canonical to a microcanonical ensemble). In the 
last step, a more convenient expression is developed by finding the linear ODE satisfied by 
the Taylor series and employing the method of dominant balance. After the main result is 
derived, some aspects of the derivation are examined further. In particular, the specialization 
to one spatial dimension deserves its own discussion. There is also a generalization of the 
Cardy formula for arbitrary dimension proposed in 2000 by Erik VerlindejOJ- Determining 
the scope of its validity has received a large amount of attention |1 7H2 1] . The extent to which 
our results agree with this Cardy- Verlinde formula will be discussed at the end. 
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II. CALCULATING THE DENSITY OF STATES 



A. The partition function 



A field theory can be regarded as a collection of degrees of freedom labelled by a momentum 
p and possibly other quantum numbers labelling the spin and/or charge. When a particular 
degree of freedom having squared momentum p 2 is excited, the energy added to the system 
is a/p 2 + m 2 by the relativistic dispersion relation. In a bosonic theory, this excitation can 
be repeated infinitely many times, so if we restrict our attention to this degree of freedom, 
the single mode partition function is: 



Z(p 2 ) = 1 + e ~l 3 \ / P 2+m2 + e -2/Vp 2 +m2 + . . . 

= (l - e-^P'+^y 1 (4) 

How many degrees of freedom have a squared momentum of p 2 l This question depends 
on the manifold, M. Allowed momenta are eigenvalues of the Laplace-Beltrami operator. 
This is a map between the two function spaces A : C 2 (M) — > C°(M) where the subscript v 
indicates Dirichlet boundary conditions. The eigenvalue equation is: 

A/ = -Xf 

where A = p 2 . Understanding the spectrum of this operator is related to Kac's problem; 
"Can you hear the shape of a drum?" Weyl's formula - one of the earliest positive answers to 
this question - will be immensely useful to us. Weyl's formula states that #(A), the number 
of Dirichlet eigenvalues up to A is given by: 

#(A) = (2ir)- d u d V\% + O (a^- 1 ) 

where V is the volume of the manifold and Ud is the volume of a unit ball in M d [22j 123] . 
The O indicates that this asymptotic formula for the number of eigenvalues becomes 

increasingly accurate when A is large. By differentiating this quantity, we can determine 
g(\) - the number of eigenvalues between A and A + dA. 

^ (A) = dA 



At this point, it is worth clearing up a potential source of confusion. While g(X) gives the 
degeneracy of the eigenvalue A and A is a way of labelling energy, g(X) is not the same 
as p(E). In single particle problems, Weyl's formula is sometimes referred to as a formula 
for the density of states[2l], but in this problem, we are interested in the density of states 
for the entire Fock space. We will therefore use the term degeneracy for g(X) and reserve 
the term density of states for p(E). While this formula for g(X) gives us a degeneracy for 
the eigenvalue A coming from geometrical considerations, we know that a field theory can 
introduce further degeneracy as part of its description. Therefore, if s is the number of 
internal states, there are really sg(X) ways of exciting a field to an eigenstate that has a 
squared momentum of A. These excitations contribute a factor of 



to the partition function. The full partition function we seek is the product of Z(X) over all 
eigenvalues A. Taking the log of a product turns it into a sum, so log Z = sg(X) log Z{X). 
In the large volume limit, the spectrum of the Laplacian becomes continuous, meaning that 
our sum over A becomes an integral. 



In the second last step, we have Taylor expanded the logarithm and in the last step, we have 
used the fact that A = p 2 . To proceed further, we must make the field theory conformal by 






setting the mass to zero. This allows us to use the identity J °° e 



cx x d 1 dx 



{d-l)\ 



\ogZ 




(27T(3) d ^ n d+l 

v ' n=l 

sd\u d £(d+l)V 

(27T/3) d 



(5) 
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The analysis so far has been applied to a bosonic theory, but very little changes when 
applying it to fermions. Each mode can have or 1 excitations so instead of (Q, we have: 

Z(\) = l + e- f)VX (6) 

Proceeding to calculate ^ in the same way, we get: 

POO 

logZ = / sg(X) log Z(X)d\ 
Jo 

/oo 
log (l + e-^Ai-MA 



sdudV 



2(27r) d 



(2vr) c 



sd\u d V ^ (-1) 



E 



n+l 



(27r J 9) d ^ n d+1 



(7) 



(27T/3) d 

a 

where C*( cr ) = (1 — 2 1 ~ CT ) £(cr) is the alternating zeta function. The functional form Z = eP d 
will be essential in the next step where we link the partition function to the density of states. 



B. From Laplace-Beltrami operators to Laplace transforms 



We are now interested in finding p(E) from Z(f3). By definition, the partition function is 
the Boltzmann factor e~ l3E summed (or integrated) over all states. Grouping the Boltzmann 
factors by energy makes the following formula apparent: 

POO 

Z= p{E)e-P E dE 
Jo 

This is recognizable as the Laplace transform of the density of states so it is natural for us 

a 

to consider the inverse Laplace transform of ef d . The following manipulations rely on a few 
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straightforward identities involving hyperbolic functions. 

I re+iT 

p(E) = lim / e^e^dp 



lim 

27TZ T-s-oo 



e+iT 



e-iT 



f3E 



1 + 



COth (^) 
i rco i pe+iT 

= — / e i/3E dp + — lim / 

27T Til T^oo J t 



I3E_ 



dp 
1 



= 5(E) + — lim 



711 T-5>oo J £ _ iT 



e+iT 



'e-iT 

tanh 



coth ( 



-dp 



26" 



1 — tanh 



(2?*) 



d/3 



= 5(E) + — lim 



e+,T tanh (^)+ tanh 2 (^) 



711 T-5>oo J e _ iT 
e+iT 



i re+11 

= 5(E) + — lim / 

7TI T-,00 J e _ iT 



1 — tanh 



sinh 



2/3 c 



2/3 d ) 

+ sinh 2 



d/3 



dp 



To proceed further, we will Taylor expand e l3E . The interesting part of the density of states 
(the part without the delta function) now splits into two pieces. 



x E k 



fc=0 



lim 



1 

2tH 



e-iT 



/3«sinh(-)d/3 



711 



+%T 
e-iT 



P k sinh 2 



2/3 c 



d/3 



(8) 



To evaluate the first integral, we will make the substitution a — ' 



lim 

2ni t^-oo 



t^Tt? sinh (aa d ) 



k+2 



da 



The integrand has a pole at ct = and the contour over which we are integrating is a vertical 
line, slightly to the right of this pole, with a vanishingly small height. 
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Im(a) 



lm(a) 



Re (a) 



Re (a) 



(a) Closed contour 



(b) Open contour 



FIG. 1: The integration to consider in the complex plane. 



The figure on the right shows the desired contour, but fortunately, the integral over the 
contour shown on the left is equal. The leftward and rightward pieces of the contour cancel 
out while the downward piece is at infinity and does not contribute. By the residue theorem, 



— - lim 

2lTi T-S-OO / e-iT 



sinh (aa d ) 



a 



k+2 



da 



( sinh (aa d ) 
Res ^^;0 



a 



k+2 



ifc-d+1 



sinh (aa d ) 



(k - d + 1)! da k ~ d+1 a d 

1 d k - d+1 a^a 2 ^ 

(k - d + 1)! da k ~ d+1 ^ (2n + 1)! 

\ / n=0 



a=0 




Odd 

a 

otherwise 



The second integral in ([8]) can be computed using the same substitution and the same 
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contour. All that changes is the parity of the residue. 



1 r-$§x sinh 2 (f a d ) 
— lim / A4 — -da 



TTl T-S>00 / e-iT 
e 2 +T 2 



a 



k+2 



sinh 2 f§a d ) 
2Res — ^; 



a 



k+2 



ifc-2d+l 



sinh 2 (f a d ) 



(k - 2d + 1)! da*- 2d+1 
1 



a 



2r/ 



o=0 



(A;-2d+l)!da fc - 2d 



d k-2d+l ~ a 2n a 2d(n-l) 
1 



71=1 



(2n) 



o=0 



We can now substitute both of these back into l\8 




fe+i 



even 



otherwise 



p(E) = 8(E) + £ ^7^+ E 



*±1 odd 



<w + E 



k\ (*±1)! 



k\ (*±*)! 



ri 



di kl (Mi), 



<w + E 



(9) 



j= -i!(d7--l)! 

When d = 1, this series can be expressed as a modified Bessel function but for other values 
of d this is surely something much less ubiquitous. 



C. An asymptotic series from a Taylor series 



In the previous step, the density of states we seek was expressed as a Taylor series, but 
not one with an elementary functional form. To better understand this Taylor series, we 
will define the function f(x) = J2m=i^"iX m where b m = m! ( dr ^_ 1 ) ! - The recurrence relation 
satisfied by the b m is m(dm — l)(dm — 2) . . . (dm — d)b m = 6 m _i. Multiplying both sides by 
x m and summing, 

oo 1 

V -(md) ...(md- d)b m x m - b m x rn+1 = 
' d 

m=l 

If x m were x dm , we could differentiate it d+ 1 times to turn it into (md) . . . (md — d)x dm ~ d ~ 1 
and then multiply by x d+1 to get it into the form above. This suggests that our differential 
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equation satisfied by / is: 



1 d d+1 
1 _a_ d d+1 



d 



X d+i . 



dxd 



d+1 



1 — xf(x) 

y=xd J v > 

f(x) - xf(x) 








(10) 



For a fixed d, this equation can be expressed in a more standard form. The first three of 
them are: 



d 



x 2 £,f(x)-xf(x) = 



dx 
„3 d 3 



./•(•'•) + <i.'- 2 ,;:,/(.'•) - xf(x) = o 

27x i £ [ f(x) + W8x 3 £f(x) + Wx 2 £f(x) - xf(x) = 



It can be seen without much trouble that when a power of x appears beside a derivative, 
the exponent is equal to the order of the derivative. While it is rather technical to find the 
asymptotic series for f(x), it can easily be converted into the asymptotic series for p(E) 
because p(E) = j^f(aE d ). 



Fact 1. The coefficient on the highest order term in (10) is d 



Proof. This is clear because ^+\ f(y d ) = j^s [dy d 1 f'(y d )] ■ Every time one differentiates /, 
one gets another factor of dy d ~ x . Doing this d more times yields d d+l y d2 ~ l f (d+1) {y d ). Letting 



xd and multiplying by kx^r , we get d d x d+1 



x 



□ 



We will now set f(x) = e 50 ^ and observe that the highest order derivative gives us 
the term d d x d+1 (^r) d+1 after we cancel e 50 ^. The other terms will all be of the form 
(^) ni • • • xi where n 1 + 2n 2 + ---+ j nj = j and j < d + 1. 



Fact 2. Choosing So(x) ex 

r d+l fdS \ d + V 



x d +! caus 



es all such (^J)™ 1 



O X 



Axi 



to be 



Proof. Substituting our guess for Sq(x) into the term above, we find that it is proportional to 

x i = x 3- ni (4r)- na (3¥r+ 1 )— - n 4dTT+^- 1 ). We need to see if this 



d 



d 

. d+1 



+J-1 
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is o(x). Clearly it is because to maximize j — n\ (^tj) —n 2 (^rj + l) rij (^rj + j — l) , 

we only want to be subtracting the n\ part, so we let all other n« = 0. We now want to 
minimize j (l — ^) and the way to do this is to let j = d + 1, giving us x. Anything less 
will give us something o(x). □ 

With the knowledge that we can make all terms but d d x d+1 (^f) d+1 negligible, our differen- 
tial equation becomes: 

d d x d+1 (^)" +1 = * 
dSo - 1 -x^-i 



dx dd+i 

So(x) = (12) 

The solution to the differential equation when we drop all but one term, is precisely the 
solution that we have shown in Fact 2 to be consistent with the dropping of the terms in 
the first place. Thus, to zeroth order, the asymptotic behaviour of f(x) must be log/(x) ~ 
( ^ d+ <i d + i ' With a little more work, we can figure out what it is to first order. 



Fact 3. The coefficient on the second highest order derivative in (10) is d 



Jd d 2 -l 



Proof. Again, starting with ^yd+\ f(y d ) = ^ [dy d ~ 1 f'(y d )] , we have to differentiate y d ~ 1 
once and differentiate / every other time in order to get an f^ d \y d ). If we wait until we have 
4- [d d y d ( d ~^ f( d \y d )~\ , the last derivative gives us d d d(d-l)y d ^ d -^- l f d \y d ). If we were to dif- 
ferentiate the factor of y earlier when it had y T ^ d ~ l \ we would get d d r(d-l)y d ^ d -^~ 1 f i - d \y d ). 
Adding these up, we get d d (d-l) fe=i r ) y^'^ 1 f [d) {y d ) = d d+1 ^y d2 - d ~ 1 f^(y d ). Using 
y = x^ and postmultiplying by hx^r again, this becomes d dd -=^x d f^ d \x). □ 



Since we have a solution for Sq(x), we will now let f(x) = e s ow+Si(x) _ when this is inserted 



into (10), the highest power of x multiplying a single derivative of S\(x) will come from the 



first term in the expansion of 



d+1 d_ dSi 



dd d+1 - - - ■ - + ] (13) 
\d + l ddi dx 1 



d+1 
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and it will be d (^i)^ 1 ( + ' ) " +1 x d+1 ~^ Also, by looking at the zeroth term 



in the expansion of (13), we see that an x exactly cancels the — x present in the definition 



of (10). 



Fact 4. Choosing S\(x) ex logx makes x d+1 d + x (^r) = £ d+1 the dominant term. 

Proof. Recall that when So(x) was the exponent, we had terms that looked like 
(^if) 1 • • • ( ) X ^ wnere n\ + 2^2 + • • • + jrij = j and j < d + 1. This time, we replace 

1 / 1 \ ^1 / 1 ■ \ n 3 

S (x) with x« +S'i(x) so this product looks like fx^n -1 + ^ J ... (x^+ I ~ j + ^Jm x J . 
Terms in the expansion of this product look like: 



dx / \ dx J 

ryr % 3 x d+1 3 -ki-2k 2 jkj - (m -hi )- 2 (n 2 -& 2 ) j{rtj~kj) 

fexH hfc,' 



x d+1 



To maximize the exponent, we must turn it into an d the way to maximize this is 

to let Hi — j — d+1. However, we have aleady mentioned that the x term should vanish so 

d 

the next highest power of x we can get is i^ 1 . We realize this either with ni = j = d or 
with j — d+ 1 and k\ — n\ — 1. Everything else is o ^^j. □ 



This allows us to drop most of the terms in (10) just like before. Specifically, we will drop 



all powers of x, lower than x d + x and find that this is consistent with having Si(x) oc logx. 

d 

The next two facts are devoted to solving for the coefficients that appear beside x^ 1 . 
Fact 5. ^e ram can be written as x~ n e rxm Y^j=o a j( r x m y ■ 

Proof. We are dealing with 4^=% (mra; m_1 e r ' x ' m ) . If we differentiate the exponent k times, 
we get m k r k x k( - m ~ 1 ^e rxm . If we then differentiate the power of x n — k times, our result is 
proportional to the exponential times r k x k ( m - 1 )-( n - k )= x ™(™ m ) fc _ □ 



In the situations where we need to apply this fact, m = -rj and r = One way 

d i i 

to get an x^ 1 term is to let n = j = d. In this case, the coefficient we find is r m = 



14 



f^^a — J d+1 (d+i) j ■ We must remember that this is multiplied by the coefficient on the d th 
derivative found in Fact 3. 



Fact 6. The other route to an x term - choosing n = d+l and j = d - yields a coefficient 
d? 1 



of 



2 J+ ■ 
d<i+i 



Proof. We have £^ ( 



mrx 



e ra ' m ) and we have to differentiate the exponential all but one 



time. Waiting until 4^ ( m d r d x d ( m ^e™), this differentiates to d(m — l)m d r d x dm d 1 e ra ' r 



Had we done this for x s ^ m ^, we would have had s(m — l)m d r d x 



d^d^dm—d—l 



times the exponential. 



We can therefore add these up to get ^4- (^zi)' 



1 \d ( {d+l) d + 1 \ d+l 



Ed 
s=l 



-d 2 1 

2 d 2 
dd+i 



□ 



We must also remember that this is multiplied by the coefficient on the (d + l) st derivative 
found in Fact 1 . The following equation is the result of keeping only the terms found to be 
significant in Fact 4 complete with prefactors. The left hand side uses the coefficient observed 



from (13) and the right hand side uses the two coefficients found above corresponding to the 



two ways of constructing . 



d 



d + l 



d-l 



{d+l) d+1 
d~ d 



d 
d+l 



14-— 
X d +! 



d^l 

dx 
dSi 
dx 

Si{x) 



1 

x 



9 d 2 
' d 2 



d 



d 2 



l l 



2 

d-l 



d 2 
d d +^ 



2(d+l) 2 

log x + log C d 



Finally substituting (12) and (14) into f{x 



2{d + l) 

e s (x)+Si(x) we obtain the result: 



1 f(^) d+1 V +1 x5+ T 

f(x) ~ C d x*V+Ve\ d 



p{E) ~ C d a 2l > d+1 )E 2( - d +Ve 



t r lj-T, d + 1 \ d+l 1 _ _ 



d 

X d+1 



(14) 



(15) 



A remaining point to address is the integration constant Cd- We used the method of dominant 
balance to find the asymptotic behaviour of /, a function defined by a Taylor series. However, 
what we really did was find the asymptotic behaviour of any function that solves the same 
differential equation as /. This equation is linear, so there is no way that such an analysis 



can tell us the constant by which the expression ( 15 ) should be multiplied. The integration 



15 



constant can be fixed if we substitute the same large number into /'s Taylor series and 
/'s asymptotic series and demand that these expressions be equal. This has been done for 
several values of d in the plot. Looking back to Q, it is no coincidence that our computed 
value of C\ is very close to 



0.30 




FIG. 2: The overall constant appearing in the asymptotic series for the density of states as 
a function of dimension. 



It is also worth noting that this method could in principle be continued to find additional 



terms in the asymptotic series for p(E). The only differential equation (10) we could write 
down for a general d had coefficients that were not very explicit. This fact would make the 
higher order terms progressively more difficult to compute. However, for a fixed d, there is 
no reason to stop at the f(x) = e s ow+Si(x) ansa ^ z p or d — 1, it is a simple exercise to find 
the more accurate asymptotic expansion 



1 2y^ 



f(x) ~ dx^e v ^ 
p(E) ~ Cia^E~^e 2y/ ^~^v^E (16) 
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which comes from considering f(x) = e s o( x )+Si(x)+s 2 {x) w ] ien Sq(x) and Si(x) are already 



known. 



D. Limitations 



Going back to our calculation of the partition function, a formula we used was: 



\ogZ 



s^doj d V 
(2n) d 




oo 



(±1)*+1 



n 



e - n p Pp d-l dp 



The approximation of a continuous momentum spectrum is only useful if the large momenta 
contribute significantly to Z. This contribution vanishes as we raise the value of /3 so our 
formula for the partition function is only valid at high temperatures. Another way to say 
this is that we want the correction to Weyl's formula to be small. Weyl's formula for #(p 2 ) is 
roughly a degree d polynomial in p. The term that we have been using is proportional to the 
volume of the manifold but there is a second term, found by Ivrii[22], that is proportional 
to the area of its boundary dV. 



This would give the log of the partition function a -^=x term in addition to the that we 
found. Making this extra term small requires that the temperature be larger than the inverse 
length scale 

This high temperature constraint on the partition function translates to a high energy con- 
straint on the density of states. It is easy to see this if we use the fact that the peak energy 
in a canonical ensemble goes as VT d+l . Another justification comes from the formula for 
the density of states. j3E must be on the order of 1 if we want the "incorrect parts" of the 
partition function to contribute negligibly to the inverse Laplace transform. If E ^ were 

d d 

not satisfied, the p(E) cx e Ed+1 result that we found would in fact be a convolution of e Ed+1 , 



e E d , e Ed 1 , etc and the error arising from not computing this convolution would dominate 
the effect of the logarithmic correction. 

One can imagine a procedure which "squishes" the manifold such that dV increases but V 



#(A) = (2n)- d u d V\l + ^Tty^co^dVX^ 1 + o (a^ 1 ) 



d-l 
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stays the same. This will raise the energy required to make sense of our results, but as long 
as BV stays finite, such an energy will still exist. A different problem occurs if we take BV 
to zero. E 3> ^ is only the right condition if the dominant correction to Weyl's formula is 
proportional to BV. For a boundaryless manifold (such as a sphere), E > is not sufficient 
to make our main result correct. When we discuss the Cardy-Verlinde formula, we will show 
that correct condition for a sphere of radius r is E ^> r _1 . 

Knowing this, it is perhaps not useful that we used the method of dominant balance a third 



time to derive (16). We can Taylor expand it to get: 



I 7-1 \ i 7-1-3 2VaE 1= 

p(E) ~ C\a±E ±e 



C ia *E--*e 2VaE [1-0 



1 



'tE, 

In deriving this, we used the fact that Er ^> 1. However, the Taylor expansion tells us that 
the only term we derived by going beyond our main result, is one that we already assumed 
to be zero. By comparison, the additional assumptions required to apply our result to an 
asymptotically free theory, will not lead to their own inconsistencies. If E ^> Aqcd, we 
still only need to make sure that Er 3> 1. The formula for p(E) will not yield any terms 
proportional to *i . 

III. DISCUSSION 

A. Other approximations 



The calculation of (15) was split into three steps but one may ask if all of these steps are 

a 

really needed. Most of the progress is made in the first step. Once the formula Z = e$ d is 
derived and the constant a is known, there are two simple procedures that can be carried 
out to get an approximate density of states. Unsurprisingly, they give the same result. 

The first procedure is to calculate the entropy in terms of temperature using 

BE B 
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One can then solve for the energy in terms of temperature using E = -^J^. This is really 
the expected energy of the canonical ensemble, but energy fluctuations away from this are 



surpressed in the thermodynamic limit. Inverting this, one can insert temperature into (17) 
to find S(E). 

~ d 

The second procedure is the one used by Cardy|7J. The power law S = CE^ 1 for some 
prefactor C is known by conformal invariance. The definition of temperature states that 
§§[e = which can be used to solve for the expected energy E in terms of (5 and C. In the 

a 

saddle point approximation, Z = e s ( E o)-PEo _ Demanding that this be equal to e& d allows 
one to solve for C and thereby solve for S(E). 

What we get from both these analyses is that: 

S = rr a d +^Ed+i 

\ d* J 

If we exponentiate this, we get the zeroth term in the asymptotic series for p(E) - our 



result (15) without the logarithmic correction. The role of the logarithmic correction can 
be understood from dimensional analysis. Boltzmann's formula states that S = \og(Q(E)) 
where Q(E) is the exact number of states that have energy E. Another way to say this is that 
Q(E) gives the number of states between E and E + dE. However, the quantity Q(E) is not 
very useful for the techniques that approximate the spectrum of a thermodynamic system 
as being continuous for large energies. What we seek instead, the density of states, gives the 
the number of states in [E,E + dE) as p(E)dE which requires it to have units of inverse 
energy. It is essentially Q(E) times an inverse energy prefactor. In some sense, the second 

1 -d-2 

and third steps in our main calculation exist to solve for this prefactor as a 2 ( d+1 > E 2( - d+1 "> . 

Even though we focused on finding the microcanonical density of states, our analysis is 
similar in spirit to two studies of how the entropy fluctuates in a canonical ensemble [251 126] . 
These authors also express the density of states as the inverse Laplace transform of the 
partition function. Rather than inserting a particular form for the partition function, they 
leave it general in order to obtain the entropy's deviation from S(f3 ) in terms of S(f3 ). The 
saddle point approximation is used to evaluate this inverse Laplace transform so that it has 
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an exponent of logZ(/3 ) — 0qE. Both papers agree that the fluctuation is: 

S = S - - log Cy 



c v = ^ Z 



d(3 2 

However, they use different methods to calculate Cy and [25] makes note of a conflict: [25] 
finds S 1 = if>(/?o) — llogS^o) while |26j finds S = S((3 ) — |logS'(/3o)- Since we have 
an explicit form of the partition function, we can see which result is reproduced in a free 
theory. The equilibrium inverse temperature we find from Z = e? d is p = (f ) d+1 . By 
differentiating logZ twice and using this temperature, we find that: 

C v = (d+ l)(<ia) 3 T T E^Ti oc E^h 

Since this is proportional to S((3 ) and not S((3 ) 3 , entropy fluctuations in this setting appear 
to be in favour of 1251. 



B. Relation to CFT 2 



The real strength of Cardy's approach lies in a duality between the low and high temperature 
limits of a CFT partition function. 

The partition function for a CFT on a torus is usually expressed in terms of the modular 
parameter t. A parallelogram is defined by two linearly independent vectors, and a torus is 
defined by a parallelogram if we identify the edges. If we express the two vectors as complex 
numbers, r is defined to be their ratio. This means that r is dimensionless. A generalized 
definition of the partition function is 

Z = Tr [ e -V(# Im(r)-iPRe(T))j 

where H and P are the energy and momentum operators respectively and V is the spatial 
volume (one of the periods of the torus). This becomes the regular partition function from 
statistical mechanics if r = y. Unitary conformal field theories possess the property of 
invariance under the modular group PSL(2, Z). Modular transformations take the form 

ar + b 
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where a, b, c and d are integers satisfying ad — be = 1. They are generated by two transfor- 
mations: 

S:t4-- T : t ^ t+ 1 
r 

To discuss modular invariance in the context of our result, we should consider what the low 
temperature limit of our partition function should be. At high values of (3, all Boltzmann 
factors other than the one coming from the ground state are exponentially surpressed. The 
partition function is therefore dominated by the energy of the vacuum which is — y-^p- 
Here, c and c are the left and right central charges. Since they are equal in the theories we 
consider, we will call this — f^- A proper partition function should therefore satisfy: 

Z ->■ as (3 ->■ oo (20) 

Clearly, the partition function we wrote down earlier does not have this limit because we 
assumed high temperatures in its construction. However, this limit is all that we need to 
know to derive the high temperature limit as well. Inverse temperature and time are related 
by a Wick rotation so a conformal field theory in one spatial dimension can be regarded as 
a Euclidean CFT on a torus. 



P 



2n 



2n 



4n 2 /p 



FIG. 3: A CFT in volume V at inverse temperature (3 and a CFT in volume V at inverse 

Yl 



temperature K- are related by a conformal transformation. 



Scaling both lengths by ^ turns a V by (3 torus into a j by V torus. This scaling is 
a conformal transformation and a Euclidean theory is agnostic to whether the first or the 
second co-ordinate carries temperature information. This tells us that we must have: 

P 

This provides the desired link between the high and low temperature partition functions: 



Virc 



Z ->■ as (3 ->■ (21) 
21 



The low temperature limit (20) holds in any number of dimensions with any type of in- 



teraction. The high temperature limit (21) holds with interactions as well but only in two 
dimensions. The reason for this is clear from the derivation - the number of dimensions with 
and without temperature information have to be the same. The figure shows the common 
idiom (5 i— > ^j-. With r = y, this is equivalent to switching f3 and V. In other words, this is 
an S transformation. In many cases, invariance of the partition function under S is simply 
called modular invariance [81 127] because the effect of the T transformation is trivial. 



Naturally the partition function given by ^ has the form of (21). Substituting d = 1, we 
have: 

Tip 
VlTS 

= ~w 

Even though our derivation required no knowledge of the Casimir effect, we confirm that 
the central charge of a bosonic field is equal to the number of internal states[27j. This also 
holds for Dirac fermions which have c = 1[28]. The number of spin states is s = 2 and 
the required factor of | comes from the zeta function identity C*(<x) — (1 ~~ 2 1 ~ cr ) C(er). It is 
therefore clear that if we choose d = 1, our formula for the partition function - and therefore 
the density of states - becomes valid for an interacting theory as well. As stated in the main 
text, the choice d — 1 also renders the differential equation exactly solvable so there is no 
need for asymptotic analysis. The Taylor series ^ becomes the Taylor series of a modified 
Bessel function giving the result of [8]. 



C. The Cardy-Verlinde formula 



When the temperature is large compared to the inverse length scale, what we found to be a 
negligible effect of topology on the partition function is not surprising. In [29], the authors 
state that this behaviour is fully expected. They also give an example appearing in the 
AdS/CFT correspondence PP which shows that the assertion may fail at low temperatures. 
As a consequence of Weyl's formula, we have seen that the leading term in the free energy is 
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indeed independent of the manifold. However, the sub-leading behaviour need not be. The 
sub-leading free energy on a manifold with a boundary will look much different from the 
sub-leading free energy on a manifold without a boundary. This is why our final discussion 
will specialize to a sphere of radius r. As we will see, specializing to a manifold with a 
boundary instead would cause the Cardy-Verlinde formula to break down. 

The form of the Cardy formula that we have mentioned multiple times is: 

S = 




3 

where r is the radius of the circle on which the CFT is defined. On a rf-sphere, the Cardy- 
Verlinde formula[6] looks very similar: 

1i\r 



S=-j=VEc(2E-E c ) (22) 

An entropy proportional to \/~E rather than would obviously be not be correct so some 
care has to be taken in how the Casimir energy Eq is defined. Verlinde arrives at the formula 
using the relations: 

E = E-& + -£?c 



E E = —S 1+1 * 



Attt 

B , 

c 



E c = —S 1 '^ 



2irr 

Ee is the extensive part of the energy, something that we have been careless enough to have 
simply called the energy all this time. Its dependence on 5 ,1+ s should be familiar by this 
point. The non-trivial claim is that E — E$ oc S 1 ^^. 

Taking S oc VT d , it suffices to show that Ec oc T . For d > 1, this contribution vanishes 
at zero temperature. It is therefore completely unrelated to the vaccum energies that are 
typically calculated using regularization techniques. Rather, in a free theory, it is related to 
the finite size effects in Weyl's formula. We have seen that the first term in Weyl's formula 
contributes VT d+l to the free energy (the extensive part), while the second term contributes 
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dVT d . This means that in order for the terms in the free energy to jump from T d+1 toT d l , 
the manifold must be boundaryless. 

The Cardy-Verlinde formula has been most successful at describing strongly coupled theories 
with AdS duals. The main work studying its use in a free theory is a 2001 paper by Kutasov 
and Larsen[T7]. They compute exact partition functions without passing to a continuum 
limit for free theories on S 1 x § 3 and S 1 x S 5 (some of their results can be checked in [5]). 
Most of their results show the desired T d ~ x scaling and this is because they use a shape with 
an empty boundary as their spatial manifold. By examining spherical harmonics, we can 
perform a calculation in the continuum limit that agrees with the results of [T7| except for 
the d = 3 scalar with conformal coupling. 

The <i-dimensional spherical harmonics obey the eigenvalue equation 

r 2 AY h _ ld = -l d (l d + d-l)Y h _ ld 

and have integer indices that satisfy < \lx\ < I2 < • • • < /d|3Q]. If the condition were 
h < ■ • • < Id, we could use the formula for /^-simplex numbers to solve for the degeneracy 
of each eigenvalue as (^^l^ 1 ) • However, li is allowed to be negative. Multiplying by 2 and 
subtracting the number of l 2 < ■ • ■ < Id choices so as not to double count l\ = 0, we find a 
degeneracy given by: 

n\ nf l d + d~l\ (l d + d-2\ 2l d + d-l fl d + d-2^ 
^ = \ d _i )-{ d _ 2 )= JZ\ { d _ 2 

Using the eigenvalue equation and the degeneracy, we may write: 

2l d + d - 1 (l d + d - 2 





- i; 


j2 _|_ ^jylipl 




d- 


- 1 


V(d- 


- r 






d- 


- 1 


2r 2 p 




+ d-2\ 


d-1 




d-2 / 



<i-i v d-2 j dld ~ j^i v d-2 r h 

« /(d- IV + drW n d + d- 2\ 2r 2 p 

dp 



d-2 J ^/(d- l) 2 + 4rV 
dp 
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Carrying out the same procedure as before, we have: 
>o g ^-^p p (^- 2 ),o g(1 - e -)d P 



2s 



/■oo 

/ r \l d + d-2)(l d + d-3)...(l d + l) log (1 - e-^) pdp 
■ Jo 



(d-l):j 

2s roc „2 



(d -1):, 

2s 



j jf ^(2/ d + 2d - 4)(2Z d + 2d - 6) . . . (2/, + 2) log (l - e"^) pdp 
(^4rV + (d - l) 2 + d - 3)( v / 4rV + (d - l) 2 + d - 5) 



(d-l)!7o 2 rf - 2 
. . . ( v /4r 2 p 2 + (d-l) 2 + 3 - d) log (1 - e-P p ) pdp 

Notice that the product in the integrand is already factored as a difference of squares. If 
d = 2, it is the empty product. If d is an even number greater than 2, it is a product of 
(4r 2 p 2 + (d - l) 2 + (d - 3) 2 ), (4rV + (d - l) 2 + (d - 5) 2 ), etc. If d is odd, there is an 
extra ^ Ar 2 p 2 + (d — l) 2 left over which we Taylor expand as 2rp+ ^ d ~^ . Collecting all the 
highest powers of r in the integrand, we get r d p d ~ x . We would now like to find the term with 
the next highest power of r. In the even case, it is: 



_ r d-2 p d-3 



d-4 

(j-l) (d-l) 2 -X>+l) 



r d - 2 p d - 3 (d 3 -Ad 2 -lid- 6) 
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In the odd case, it is: 

\ r d~2 p d~3 



d-3 

(^-l) (d-lf-±(2k) 
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r d - 2 p d - 3 d(d 2 - 3d +11) 



We can therefore evaluate 



9s f°° 1 
logZ = -T^yy jf (rV- 1 + -Md)r rf -y- 3 ) log (1-e-^) dp 



2dsr d ((d + 1) sh(d)r d - 2 ((d - 1) 



+ 



6(d- l)(d-2)/3 



d-2 



where 



/i(d) = < 



d = 2 

d 3 -Ad 2 - lid- 6 d even 
d 3 -2d 2 + lid d odd 



(23) 



25 



Again the zeta function is the only thing that would change if we had done the calculation 
for fermions. Notice that 6 is a root of h(d). We now have the extensive and sub-extensive 
parts of the free energy and the entropy: 



drrd+l h(d)[s((d - 1) + s*C(d - 1)] 



F = -2d[<(d + 1) + s*C(d + l)]r d T d+1 - 6 (d-l) d- 2) 



S = 2d(d+l)[sad+l) + s*C(d + l)}r d T d + h{d)[sad Jjj^A^ 1)^ -2^-2 

Defining the Casimir energy as E + dF (the deviation from the Euler relation [6]), it is 
straightforward to compute the energies: 

E E = 2dM(d + 1) + sT(d + l)]r*T" + " gggjj^ " 

h(d)[sC(d-l) + s*C(d-l)] d _ A _ x 
c 3(d-l)(d-2) 

Strictly speaking, neither is a pure power of S. However, if we only require the leading terms 
to be a power of the leading term of S, the coefficients we seek are: 

A d 



47rr r (2d[s((d + 1) + s*C*(d + (d + l) 1+ 5 
B h(d)\s((d- 1) +s*C*(d- 1)1 . „ , , r>/1 . 



2vrr 3(d-l)(d-2)r 



Substituting these into (22), the corresponding Cardy-Verlinde formula is: 



3(d-l)(d-2) a C(d+l) + a *C*(rf+l) / F r9F F ; 

1 

' (2d[<(d+l) + s*C(d+l)]) 3TT r^i^ +1 



d d 

The holographic theories studied by Verlinde have = d 2 . When Kutasov and Larsen say 
that the Cardy-Verlinde formula does not hold in a free theory, they are talking about this 
more restrictive type of Cardy-Verlinde formula. The formula above holds on any sphere 
that is not a 2-sphere or a 6-sphere and as noted in [T7], the overall coefficient depends on 
the matter content of the theory. Extending this formula to a general boundaryless manifold 
would require us to know the third term in Weyl's formula. 
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